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Resumen 


BAUTISTA-DE-LOS-SANTOS, Q.M. Determinación de 
caudales ambientales en la cuenca del río Yuna, República 
Dominicana. Tecnología y Ciencias del Agua. Vol. V, núm. 3, 
mayo-junio de 2014, pp. 33-40. 


Se han estimado los caudales ambientales utilizando el 
método del análisis del Rango de Variabilidad (RVA, 
por sus siglas en inglés) y los Indicadores de Alteración 
Hidrológica (IHA), en tres puntos de la cuenca del río 
Yuna (Embalses de Hatillo, Rincón y Blanco), localizada 
en República Dominicana. Las series de caudales diarios 
para cada sitio de estudio fueron obtenidas mediante una 
combinación de generación sintética de caudales mensuales 
y desagregación en caudales diarios. Los treinta y tres (33) 
parámetros calculados en cada punto fueron utilizados para 
prescribir regímenes de caudales ambientales, demostrando 
que la aplicación del método es factible en el contexto local, 
pese a las dificultades encontradas en el levantamiento y 
procesamiento de los datos. 


Palabras clave: caudales ambientales, caudales ecológicos, 
alteración hidrológica, Enfoque del Rango de Variabilidad, 
Yuna. 


Introducción 


Los “caudales ambientales” o “caudales ecoló- 
gicos” son regímenes de flujo capaces de 
mantener el funcionamiento, composición y 
estructura del ecosistema fluvial del cauce 
por el cual circulan en condiciones naturales 
(García-de-Jalón y González-de-Tánago, 1998). 
Una revisión global sobre el tema indicó la 
existencia de 207 metodologías de estimación de 
caudales ecológicos propuestas desde la década 
de 1940 (Tharme, 2003); estas metodologías han 
sido clasificadas y descritas por varios autores 
(Arthington y Zalucki, 1998; Tharme, 2003; 
Acreman y Dunbar, 2004). En la actualidad se 
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favorece la aplicación de métodos que consideren 
la variabilidad natural de los caudales y su 
impacto sobre los ecosistemas, en sustitución 
de los métodos inicialmente propuestos, que 
asignan un valor fijo de caudal ambiental durante 
todo el año (por ejemplo, el Método Montana; 
Tennant, 1976, citado por Tharme, 2003), el 
cual prescribe el caudal ecológico como el 10% 
del caudal medio anual y es todavía el método 
utilizado en la República Dominicana). Es por 
esto que la presente investigación se propone 
realizar una estimación de caudales ambientales 
aplicada a la República Dominicana, utilizando 
métodos más abarcadores y acordes con las 


recomendaciones actuales. Dicha estimación 
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se realizará en tres sitios de la cuenca del río 
Yuna (5 630 km?), la mayor cuenca del país, 
correspondientes a las tres obras de regulación 
de mayor importancia e impacto en la cuenca: 
los embalses de Hatillo, Rincón y Blanco, cuyas 
capacidades de almacenamiento son 441, 75.50 
y 0.73 Mm”, respectivamente. 


Metodología 


La estimación de caudales ecológicos se ha 
hecho aplicando el método del Enfoque del 
Rango de Variabilidad (RVA, por sus siglas 
en inglés) (Richter et al., 1997), la cual es 
una metodología fundamentada en el uso 
de Índices de Alteración Hidrológica (IHA, 
siglas en inglés) (Richter et al., 1996). Los 
IHA son 33 parámetros relevantes desde el 
punto de vista ecológico, divididos en cinco 
grupos de acuerdo con cinco características 
fundamentales de los regímenes hidrológicos: 
magnitud, sincronización, frecuencia, duración 
y tasa de cambio, los cuales son utilizados para 
prescribir caudales ecológicos. La metodología 
RVA ha sido ampliamente utilizada (Richter 
et al., 1998; Bragg et al., 1999; Black et al., 
2002; Brunke, 2002; Baeza-Sanz y García-del- 
Jalón, 2004, Opperman, 2006; Limbrumner, 
2010) debido a la facilidad de su aplicación y 
flexibilidad, comparado con los métodos de 
simulación de hábitats. Para la aplicación del 
método se han seguido los siguientes pasos 
indicados por Richter et al. (1997): 


1. Generación de series de caudales diarios. 
series 
de caudales diarios encontradas en las 


Para los sitios estudiados, las 
instituciones que registran este tipo de 
datos estaban incompletas y no tenían la 
extensión requerida (20 años o más), por 
lo que generaron de manera sintética tres 
serie de 50 años de caudales diarios (una 
para cada uno de los sitios de estudio) 
utilizando la siguiente metodología: 


a) Análisis de datos (estimación de las 
estadísticas básicas). 


b) Normalización de series de caudales 
mensuales utilizando transformaciones; 
las transformaciones usadas fueron la 
logarítmica, gamma, Box Cox y Potencial, 
que fueron evaluadas utilizando la Prueba 
de Filliben (los coeficientes de correlación 
deben ser mayores que 0.987 para n > 100; 
Rodríguez, 2012) y la Prueba del Sesgo 
(valores entre +0.5 y -0.5 para aceptar la 
normalidad de los datos; Brown, 2012). 

c) Ajuste de modelo ARMA(p,q) periódico 
univariado o PARMA(p,q) univariado a las 
series de caudales, utilizando el método de 
los mínimos cuadrados. 

d) Generación estocástica de serie de caudales 

PARMA(p,q) 

Fueron generadas 


mensuales utilizando el 


ajustado. series de 

cincuenta años de extensión utilizando el 

modelo ajustado para cada uno de los sitios 
de estudio; el modelo PARMA (Periodic 

Autoregressive-Moving-Average) ha sido 

ampliamente utilizado en la generación 

de caudales estacionales y caudales 
periódicos (mensuales), especialmente los 
modelos con orden bajo (PARMA(L1,0) y 
PARMA(1,1); (Salas et al., 1997). 

(e) Desagregación de caudales mensuales a 
caudales diarios utilizando la metodología 
de Acharya y Ryu (2012), la cual preserva 
el balance de masa y las características 


estadísticas de la serie histórica. 


La generación de caudales mensuales 
descrita anteriormente fue realizada utilizando 
el software Stochastic Analysis, Modeling, and 
Simulation (SAMS), versión 2007, desarrollado 
por Sveinsson et al., en el año 2007. La 
desagregación de caudales mensuales a 
caudales diarios fue realizada utilizando una 
hoja de cálculo y mediante la elaboración de 
macros en Visual Basic for Applications (VBA) 


para facilitar los cálculos hechos. 


2. Caracterización del rango natural de 


variación de caudales utilizando 33 
parámetros hidrológicos, conocidos como 


Indicadores de Alteración Hidrológica 
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(IHA, por sus siglas en inglés; Richter ef 
al., 1996). Para la estimación de los IHA fue 
utilizado un software desarrollado por The 
Nature Conservancy para tales fines (The 
Nature Conservancy, 2011). 

3. Selección de 33 rangos de variación, uno 
para cada IHA, que serán los objetivos 
del sistema de manejo. El rango de 
variación seleccionado para el análisis 
no paramétrico realizado (estadísticas 

en función de percentiles) fue el rango 

comprendido entre el percentil 25% y 75%, 

según las recomendaciones de Armstrong 

y Parker (2003). 


Resultados 


Los resultados del Análisis IHA para la presa 
de Hatillo arrojan un caudal medio mensual 
para el periodo 1950-1999, normalizado por el 
área, según se muestra en la figura 1. Se verifica 
un aumento en los caudales medios mensuales 
en los meses de mayo, y de octubre a diciembre. 
El periodo de estiaje se extiende desde enero 
hasta abril, y se verifica un periodo de recesión 
de caudales entre junio y julio. El rango de 
variabilidad de los caudales medios mensuales, 
estadísticas de caudales bajos anuales y otras 
mediciones adicionales que se utilizan en el 
RVA como objetivos de manejo se presentan 
en el cuadro 1. Para el embalse de Rincón, el 
caudal medio mensual para el periodo 1950- 
1999, normalizado por el área, se muestra en la 
figura 2. Se verifica un aumento en los caudales 
medios mensuales en los meses de mayo, y 
de octubre a enero. El periodo de estiaje se 
extiende desde febrero hasta abril, y se verifica 
un periodo de recesión de caudales entre junio 
y septiembre. El rango de variabilidad de los 
caudales medios mensuales, estadísticas de 
caudales bajos anuales, y otras mediciones 
adicionales que se utilizan en el RVA como 
objetivos de manejo se presentan en el cuadro 1. 
Los resultados del análisis IHA para el Embalse 
de Blanco arrojan un caudal medio mensual 
para el periodo 1950-1999, normalizado por el 
área, según se muestra en la figura 3. Se verifica 


un aumento en los caudales medios mensuales 
en los meses de mayo a julio, y de septiembre 
a noviembre. El periodo de estiaje se extiende 
desde diciembre hasta abril, y se verifica un 
periodo de recesión de caudales en los meses 
de julio y agosto. El rango de variabilidad de 
los caudales medios mensuales, estadísticas 
de caudales bajos anuales, y otras mediciones 
adicionales que se utilizan en el RVA como 
objetivos de manejo se presentan en el cuadro 
1. 


Discusión 
Fase de recolección/adquisición de datos 


Las series de caudales diarios disponibles 
del río Yuna tienen datos faltantes desde días 
hasta años, y datos en periodos de tiempo no 
coincidentes, lo que implica un mayor trabajo 
a la hora de adecuar dichas series para su 
posterior uso. Los archivos con las series de 
caudales son entregados en formato .txt, con 
una configuración no fija (numerosos saltos y 
espacios, en un formato no estandarizado), lo 
que dificulta trabajar con los mismos debido a 
que primero deben ser organizados. 


Fase de generación de caudales diarios 


En el proceso de generación de caudales se 
dificulta la aplicación de técnicas sencillas 
de extensión de caudales diarios, como la re- 
gresión lineal, debido a que las estaciones hi- 
drométricas tienen numerosos datos faltantes y 
datos en periodos de tiempo no coincidentes. 
La segunda mejor opción es la generación 
sintética de los caudales, técnica utilizada en el 
presente estudio para la obtención de las series 
de caudales diarios, pero que tiene limitaciones 
en la generación de caudales bajos, por lo 
cual, los resultados obtenidos en la estimación 
de los indicadores, tasa de ascenso, tasa de 
bajada y número de reversas, deben tomarse 
con precaución (Richter et al., 1997). La tercera 
técnica recomendada, que es la modelación 
hidrológica, requiere de numerosos datos 
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Cuadro 1. Rango de variabilidad de caudales, embalses de Hatillo, Rincón y Blanco. Periodo del análisis: 1950-1999 (50 años). 


Hatillo Rincón Blanco 
10% | 25% | 50% | 75% | 90% | 10% | 25% | 50% | 75% | 90% | 10% | 25% | 50% | 75% | 90% 

Parámetros 

Grupo 1 

Marzo 3.829| 11.75] 19.84] 29.85| 35.07| 194] 3.875 6.731 9.635) 1448| 2.035| 2.237| 2458] 3.388| 5.539 
Abril 9.539) 14.43] 22.32| 32.25| 36.23] 2.84 5.3) 8.183] 10.88] 1281] 2.012) 2.132| 2574] 3.353| 5.839 
Mayo 10.11, 15.87| 27.2| 51:16| 80.91] 5.2821 7.178] 9.58| 16.781 229| 2311] 3.114| 4.071 7.57| 9.743 
Junio 6.725| 13.93] 19.04| 29.31] 925| 2:912| 5.978| 8.7551 13.43) 18.02| 4.483| 4716]  5.29| 7.362| 10.56 
Julio 6.065| 9.567| 18.031 29.59| 67.31 3.175 5.91 7.535 8.858] 1243|  3.559| 3.964] 5.052| 5.927| 7.463 
Agosto 9.666| 19.51] 21.82| 42.52| 67.67| 3.348] 4.343| 6.495 9.825| 1161] 2309] 2925]  3.76| 4.573| 5.331 
Septiembre 6.264| 11.55] 25.54 42.39| 65.59| 4314] 6.061 7.768, 8:984| 1046| 3.6931 4252| 5.668] 6.381| 6.921 
Octubre 15.56| 22.75] 29.711 41.84| 58.02| 4.545] 6.2905| 8.49| 11.05] 12.67| 4492] 6.045] 7.0611 8.285| 9.99 
Noviembre 7.067| 24.6] 43.831 61.03| 77.5| 3.37| 6.778 9.0551 11.28| 1749| 6.187| 6.74] 7.644] 10.03] 12.43 
Diciembre 10.83| 18.26| 38.16] 70.571 so] 186) 4405|  14| 17.43] 2154| 2.508] 2.999| 5.656 7.204| 8.131 
Enero 8.749| 13.95) 25.58| 31.81| 37.65| 2.143] 3.155 1171] 17.44] 2344|  2.37| 3.652) 4928] 6.202| 7.284 
Febrero 4.285] 12.46| 18.23| 22.86| 28.28 0.9) 2414| 7.711 10.94| 17.59| 2.534| 3.002| 4066| 4.818| 5.373 
Parámetros 

Grupo 2 

1-día mínimo | 2544| 2.798| 4291| 6.38 7806] 0.418] 0.87| 1.615| 2.113] 3.607| 1064] 141] 1673 1962| 2.358 
3-días mínimo | 2.653| 2.945] 4.406] 7.211| 10.48| 0.494| 0.8817| 1.688] 2.317| 3.976| 1.064) 1457| 17551 2.006| 2.489 
7-días mínimo | 2.72| 3.143| 4.514| 8.181| 13.21| 0.7216] 0.935] 1.834] 2.595] 4.656] 1.064) 1.484] 1.792| 2:122| 2.632 
E 3.898] 4.621] 9.249| 13.42| 19.9| 0.9857| 1.88| 2.3731 4.718| 6.422| 1835| 211] 2293] 2738| 3.335 
pa 7.024| 9.892| 16.39) 22.3| 27.71| 2.503] 3.74| 5.683| 7.912| 9.383] 2.892| 3.364] 4.079] 4347| 4765 
1-día máximo 125| 19 2418| 458.1| 852,5 27.24] 37.09| 60.29 
3-días máximo | 109.2| 142.3| 181.3] 247.9| 367.6 19.511 31.53| 36.22 
7-días máximo | 72.04) 106.11 127.9| 179.3| 230.9 14.44| 19.49| 21.77 
a 37.3 56.331 82.78| 116.8| 171.7 1133| 13.61] 15.85 
maximo 

90-días 

E 27.82) 41.66] 58.07| 78.48| 116.1 8.419) 10.73| 1211 
maximo 

Núm. días 0 0 0 0 0 0 0 0 
cero caudal 
o caudal | 0. 07558|0.09556| 0.1478| 0.2172|0.3273 0.1853 0.3732| 0.1749 0.3189| 0.381| 0.4225 
Parámetros 

Grupo 3 

Pee de 330.9 32| 84.5| 107| 2158 91 122| 123 
caudal mín 

a 160.9)  244| 305.5] 350| 26.2 270 ED) 323 
caudal máx 

Parámetros 

Grupo 4 

Conteo jptlsos 3 5 10 14| 18 9 1) 158 
bajos 

ade 2 3| 475 7 118 4 6.5 9 
pulsos bajo 

Eonico puso 6.1 9) 14  18| 227 14 17 19 
altos 

Duración 1 A 1 a 2 3 4 1.1 2 3| 3.625 5 
pulsos alto 

Parámetros 

Grupo 5 

Tasa de subida | 2.253] 3.304| 4.268] 5.881| 7.986| 0.546| 0.6475] 0.86] 1.171] 1.684| 0.4159| 0.4852| 0.5897| 0.7034| 0.8057 
Tasa de bajada | -4.325| -3.172| -1.949| -1.584| 1.118] — -0.89| -0.7387| -0.6| -0.495| -0.4005| -0.3685|-0.3366| -0.2985| -0.2413| -0.2132 
aero da 118.2)  127| 131.5| 139.31 147.9| 132.1 137| 143| 151| 1568 1121 120 126 1311 135 
reversas 
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Figura 3. Distribución de caudales medios mensuales normalizados por área. 


Sitio de estudio: presa de Blanco (1950-1999) (leyenda en percentiles). 


Fase de generación de caudales ecológicos 


El régimen de caudales ecológicos estimado 


con el método es sin lugar a dudas 
significativamente mejor que el resultado 
obtenido con el método que se utiliza en la 
actualidad en el país (el método de Montana 
de Tennant (1976), citado por Tharme, 2003). 
Una vez fijada una metodología (como lo que 
se ha hecho en el presente estudio), puede 
utilizarse para todos los estudios posteriores 
de estimaciones de caudales ecológicos con 


este método. 
Operación de los embalses 


Según los datos de operación suministrados 
por la Empresa de Generación Hidroeléctrica 
(EGEHID), los embalses de Hatillo y Rincón 
no utilizan sus desagúes de fondo, salvo en 
contadas excepciones durante el año. Por lo 
que, en el caso de Hatillo, el caudal ecológico 
manejado está asociado con los volúmenes 
turbinados, que son posteriormente vertidos 


al río. En el caso de Rincón, al existir un 
contraembalse que deriva las aguas turbi- 
nadas para el riego, el caudal manejado con 
fines ecológicos del sistema completo (embal- 
se+contraembalse) está determinado por las 
salidas del contraembalse. La hidroeléctrica 
de Blanco, según los datos suministrados 
por EGEHID, libera un caudal casi igual al 
caudal de entrada, por lo que la alteración 
hidrológica que genera sería nula. Esto se debe 
a que el embalse de Blanco está completamente 
sedimentado, por lo que ha perdido su capa- 
cidad de almacenamiento (Jiménez-Ramírez y 
Farías, 2012). 


Conclusiones 


La aplicación de la metodología RVA y los 
IHA para la determinación de caudales 
ecológicos es una opción factible en el contexto 
dominicano, ante la situación actual de 
conocimiento y aplicación limitados sobre el 
tema en el país. La misma es una metodología 


más completa y abarcadora que el método 
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tradicionalmente utilizado, y constituye una 
opción recomendada para hacer estimaciones 
de caudales ecológicos, pese a las limitaciones 
y dificultades existentes. Los resultados obte- 
nidos en este estudio pueden ser utilizados 
para elaborar un programa de operación de 
los embalses consensuado entre las institucio- 
nes competentes, para poder armonizar los 
diferentes usos de los embalses, con la nece- 
sidad de proteger el hábitat aguas abajo de 
las obras. Este programa de operación de los 
embalses debe acompañarse de un programa 
de monitoreo ecológico, para verificar la 
efectividad del régimen de caudales adoptados 
y modificarlo en dado caso que sea necesario. 
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